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SUMMARY 

An analysis method and a design method are developed for one- 
dimensional, compressible flow with friction, heat transfer, and area 
change in vaneless diffusers with arbitrary profiles in the axial-radial 
plane. The effects of mixing losses due to nonuniform flow conditions 
at the impeller discharge are not considered. In the analysis method 
the variations in fluid properties, including the velocity and flow 
direction, are determined as a function of radius for a prescribed vari- 
ation in diffuser wall spacing with radius. In the design method the 
variations in effective diffuser wall spacing and in the fluid proper- 
ties are determined as a function of radius for an arbitrary prescribed 
variation in one fluid property. For efficient diffuser designs the 
fluid property selected and the manner in which its variation is pre- 
scribed will depend on viscous flow effects that are considered in 
boundary- layer studies but are not investigated in this report. 

As a result of numerical examples it is concluded that: (l) Even 

with relatively low friction coefficients and neglecting mixing losses 
near the impeller tip , the friction losses in most vaneless diffuser 
designs are considerable, as indicated by computed diffuser efficiencies 
in the low 80 ’s, and these losses result from the usually large ratios 
of wetted surface to flow area in vaneless diffusers. (2) Vaneless 
diffuser efficiencies can be improved by increased compressor flow rates 
for a given impeller tip radius so that the diffuser walls can be spaced 
farther apart (thus, reducing the ratio of wetted surface to flow area) 
without Increasing the length of the flow path in the diffuser. 


INTRODUCTION 

In radial- and mixed-flow centrifugal compressors the vaneless 
diffuser is an annular duct (fig. 1) immediately following the impeller 
and of increasing radius in the direction of flow. The high tangential 
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velocity of the fluid entering the vaneless diffuser from the impeller 
decreases vith increasing radius and^ because the tangential velocity is 
generally the largest velocity component at the impeller discharge^ the 
vaneless diffuser is an effective means of diffusing the fluid, that is, 
of converting the velocity head to static pressure. The principle by 
which this conversion is effected is demonstrated by the case for fric- 
tionless flow in the absence of heat transfer. For this case, and 
assuming that flow conditions are uniform in the tangential direction, 
the moment of momentum of the fluid is constant so that 


q^r = constant 

from which as the radius r increases the tangential velocity 
decreases and therefore the pressure rises (assuming relatively small 
changes in other components of velocity) . 

Among the advantages of the vaneless diffuser is the fact that 
choke flow occurs only if the meridional velocity qjjj (velocity com- 
ponent normal to the annulus area) is sonic. This condition usually 
corresponds to such high flow rates that choke flow occurs in the 
impeller, instead of the diffuser as is the usual case for vaned dif- 
fusers. The compressor operating range is therefore wider with vaneless 
diffusers. 

Another, and perhaps the most Important, advantage of the vaneless 
diffuser is the fact that if the tangential velocity at the impeller 
discharge is supersonic the tangential velocity decelerates from super- 
sonic to subsonic velocities without shock losses. 

Opposed to these several advantages of the vaneless diffuser is the 
disadvantage, for aircraft propulsion, of a large frontal area. This 
disadvantage may be circumvented to some extent by the use of semivane- 
less diffusers (fig. 2) in which, to diffuse the fluid more rapidly and 
thus decrease the frontal area of the compressor, vanes are placed in 
the diffuser following a vaneless section in which the velocity is 
reduced from supersonic to subsonic magnitudes. Thus, shock losses are 
avoided by diffusing the flow to subsonic velocities in the vaneless 
diffuser and the frontal area of the compressor is somewhat reduced by 
the more rapid diffusion in the vaned section. 

In order to analyze the performance of vaneless and semivaneless 
diffusers and in order to design these diffusers for optimum performance 
(including the proper setting of the vane angles in semi vaneless dif- 
fusers), it is necessary to have adequate theoretical methods of pre- 
dicting the variation in flow characteristics through the diffusers. 
These methods should include the effects of diffuser geometry, com- 
pressibility, heat transfer, friction, and mixing losses caused by the 
nonuniform flow conditions at the impeller discharge. 
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Published work on the analysis and design of vaneless and semivane- 
less diffusers is not extensive (references 1 and 2, for example). In 
reference 1 a one- dimensional method of analysis is developed for incom- 
pressible flow with friction but no mixing losses in vaneless diffusers 
with constant wall spacing and pure radial flow. In reference 2 a one- 
dimensional method of design is developed for compressible flow with 
friction. The method assumes the flow path is a logarithmic spiral and 
neglects heat transfer and mixing losses. 

In the present report methods of analysis and design, carried out 
at the NACA Lewis laboratory, are developed for one- dimensional, com- 
pressible flow with friction, heat transfer, and arbitrary variation in 
passage height in vaneless diffusers with arbitrary curvature in the 
meridional (axial- radial) plane. The effect of mixing losses is not 
considered. In the analysis method and in general for the design method 
the flow direction, or flow path, is not specified but is a dependent 
variable determined by the solution. In the design method the variation 
in diffuser wall spacing with radius is determined for a prescribed 
variation in one fluid property. For efficient diffuser designs the 
selection of the one fluid property and its optimum prescribed variation 
will depend on viscous flow effects that are considered in boundary- 
layer studies but will not be Investigated in this report. The methods 
are an extension of the work in reference 3 for one-dimensional gas flow 
in ducts with prescribed flow direction. 


THEORY OF METHOD 

Differential equations are developed that relate the change in 
dependent variables with radius, to the design and operating character- 
istics of the vaneless diffuser. The application of these differential 
equations to the analysis of flow in vaneless diffusers and to the 
design of vaneless diffusers for prescribed distributions of flow condi- 
tions with radius is described in a later section. 


Preliminary Considerations 

Coordinate system . - The coordinate system for a point on the mean 
surface of revolution generated about the axis of the compressor by the 
center line between the front and rear shroud of the vaneless diffuser 
is shown in figures 3 and 4. The cylindrical coordinates r, 0, and 
z give the radial, tangential, and axial positions of the point, respec- 
tively. The effective diffuser height h (fig. 3) measured across the 
passage in the direction normal to the mean surface of revolution is a 
function of r only 
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h = h(r) (1) 

(All symbols are defined in appendix A.) The effective height of the 
diffuser at each point on the mean surface of revolution is equal to the 
geometric height of the diffuser minus the assumed displacement thick- 
ness of the boundary layer on the diffuser walls. Use of this effective 
height rather than the geometric height is required by continuity con- 
siderations in order to give the proper average value of the velocity 
component normal to the cross-sectional flow area of the vaneless dif- 
fuser. Only the effective height of the diffuser is considered in this 
report; no investigation is made of the boundary- layer displacement 
thickness, which can be assumed or estimated from boundary- layer theory. 

The slope of the center line between the front and rear shroud of 
the vaneless diffuser detenaines the angle a (fig. 3), which is a 
function of r only, 


( 2 ) 


Assumptions . - The principal assumptions of the analysis and design 
methods are that flow conditions are uniform across the vaneless dif- 
fuser along the height h and that flow conditions are uniform in the 
tangential direction 0. Thus, the flow becomes one -dimensional, being 
a function only of the radius along the mean surface of revolution. If 
the boundary -layer profile is ignored, the accuracy of the assumption 
that flow conditions may be considered uniform across the vaneless dif- 
fuser in the direction of h depends on: (l) the angle a, ( 2 ) the 
derivative of a with respect to r, (3) the derivative of h with 
respect to r, and (4) the ratio h/r. For values of a approximately 
equal to 90^ the assumption is accurate provided dh/dr and da/dr are 
small. For values of a less than 90^ the inaccuracy of the assumption 
will depend on the ratio h/r and the derivative doc/dr; for the 
limiting case in which h/r and dcx/dr approach zero the assumption is 
good for all values of a. In practice the values of h/r for vaneless 
diffusers are usually small and the mean shroud curvature da/dr should 
be small to avoid boundary- layer separation. Thus the assumption of 
uniform flow conditions across the passage along h should be accurate 
for all values of a encountered, except for variations due to the 
boundary- layer profile. 

The motion on the mean surface of revolution is assumed to be 
steady and, because flow conditions are assumed to be unifonn in the 
tangential direction, mixing losses resulting from nonuniform flow con- 
ditions in the tangential direction at the impeller discharge are neg- 
lected. These losses are relatively high, but experiments (reference 4, 


. -1 dr /X 

a = tan — = a(r) 

dz ^ ' 
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for example) indicate tliat they take place in the immediate vicinity of 
the impeller discharge and may he neglected in the remainder of the 
vaneless diffuser. The effects of these mixing losses can he accounted 
for (approximately) hy adjustments in the flow conditions (pressure, 
density, velocity, and flow direction, for example) at the diffuser 
inlet. 

Velocity components . - The velocity q at a point on the mean sur- 
face of revolution is tangent to the surface and has components q^, 
qg, and lii the r-, 6-, and z-dlrectlons, respectively. In 

this analysis it is convenient to consider the meridional velocity q^^ 
(instead of q^, and qg), which is tangent to the center line between 
the front and rear shroud of the diffuser in the meridional plane 
(fig. 3) and is related to q^. and q^ hy 



The flow direction 3 on the mean surface of revolution is related to 
qm and q^ hy (fig. 5) 


tan 3 = — 




(4a) 


from which 


= q sin 3 

(^4h) 

= q cos 3 

(4c) 


Fluid particle . - A fluid particle on the mean surface of revolu-- 
tion is shown in figure 4. This particle has the dimensions r dG and 
dr/sin a on the surface of revolution and the height h normal to the 
surface . 

Outline of method . - The state of the fluid at any point (r) on the 
mean surface of revolution is described by three thermodynamic proper- 
ties^ by the fluid velocity^ and by the flow direction. These five 
properties can be determined from five fundamental relations: (l) con- 

(2) ocjuilibrium in the direction of q (meridional equili- 
hrium), (3) equillhrium in the direction of q^”^ (tangential equlli- 
hrlum), (4) equation of state, and (s) the heat-transfer equation. In 
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addition to these five fundamental relations certain definitions are 
required to express the resulting equations in terms of the desired 
properties. The properties that will he used in this analysis to 
describe the state of the fluid will he the static pressure p, the 
static density p, the total temperature the local Mach numher M, 

and the flow direction p. 


Mach numher. - The local Mach numher 


M is defined hy 


M^ = 


ygE^J^T 


( 5 ) 


where 7 is the ratio of specific heats, g is the gravitational 
acceleration, is the gas constant, and T is the local, static 

temperature. From equation (5) 


M^ dr q^ dr T dr 


(5a) 


Total temperature . 


The total temperature T is defined by 

T/ 


T . = T + 
t 2Jgc 


or 



( 6 ) 


where J is the mechanical equivalent of heat and Cp is the specific 
heat at constant pressure. From equation (6) 


dr 


1 

T 


dT / ^ ^ 


1 dM^ 
W dr” 


(6a) 


From equations (5a) and (6a) 
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1 dq^ 
^ dr 


f— ) 

V * ^ *7 


1 dM^ 1 dT^ 

dr T dr 

U 


(6b) 


Fundamental Eelatlons 

Continuity . - The continuity equation for one-dimensional com- 
pressible flow in vaneless diffusers Is 


pq^jj^rh = constant 


from which 


1 dp 1 dq„ 1 dh 1 

- — + + + - 

p dr < 1 ^ h dr r 


( 7 ) 


where changes in r are understood to occur along the mean surface of 
revolution. 

Meridional equilibrium . - The equation for meridional equilibrium 
oi" ^ fluid particle (fig. 6) in the direction of q on the mean sur- 
face of revolution is obtained from a balance of the pressure forces, 
shear forces, and inertia forces (appendix B) 

, 2 „ 2 
g dp c^q cos p qg 

p dr h sin a r 

where c^ is the skin- friction coefficient. 


'Im 


dr 


■m 


(8) 


Tangential equilibrium . - The equation for tangential equilibrium 
of the fluid particle in flgirre 6 is obtained from a balance of the 
shear forces and the inertia forces (appendix B) 


Cfq P ‘i'le q^q^ 

— — = qia — + 2 

h sin a dr r 


Equation of state . - By definition a perfect gas satisfies the 
equation of state 


( 9 ) 
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p = pE*T 


(lOa) 


from which. 


1 ^ _ 1 ^ 1 dT 

pdr’~pdr'^Tdr 


(lOb) 


Heat-transfer equation . - The heat-transfer rate to the diffuser 
casing must equal the heat-transfer rate from the fluid. The heat- 
transfer rate to the diffuser casing is given approximately hy 

(T„ - Tt) ^ (lla) 


where h* is the coefficient of heat transfer, T^ is the wall, or 
diffuser casing, temperature, and dQ is the heat-transfer rate. Equa- 
tion (lla) assumes that the recovery factor at the wall is 1.0 (refer- 
ence 3, p. A-328), 

The heat-transfer rate from the fluid is given approximately hy 


dQ 


pq^ 2«rhCp 



dr 


Finally, from equations 


(lla) and (llh) 


(11b) 


Tt ar 


2h* 

pqm^iCp sin a 



(11c) 


Equation (lie) gives the change in total temperature with radius as a 
function of the heat-transfer coefficient h* . 

Reynolds * analogy . - An approximate value for h' in equation (lie] 
can he obtained from the Reynolds’ analogy between friction and heat 
transfer (reference 3) 


h’ °f 


Cppq 2 
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From which, equation (lie) becomes 

1 sec ^ /T.^ 

T.J. dr h sin a \T.j. 

Equation (12) gives the change in total temperature with radius as a 
function of the skin-friction coefficient and the ratio of wall tempera- 
ture to total temperature of the fluid. 

Beview . - A review of the theory up to this point indicates nine 
unknowns and nine equations for the analysis method. The unknowns are; 
p, p, T, T+, q^ q^j^, q^^ and p. (For the design method h is 
unknown and replaces one of these nine quantities, which is then speci- 
fied as a function of r. The angle a is a known function of r for 
both analysis and design methods.) The nine equations are: 




Equation 

Tangential velocity 

(4b) 

Meridional velocity 

(4c) 

Mach number (definition) 

(5a) 

Total temperature (definition) 

(6a) 

Continuity 

(7) 

Meridional equilibrium 

(8) 

Tangential equilibrium 

(9) 

Equation of state 

(10b) 

Heat-transfer equation 

(11c) or (12) 


The solution for the analysis method consists in combining the nine 
equations to obtain three differential equations involving three 
unknowns: T^, M, and p. These three differential equations, in 

turn, can be combined to solve, by numerical methods, for T.^., M, and 
p successively. (For the design method an auxiliary equation is devel- 

oped for i: iiS in terms of the prescribed fluid property as a function 
h dr 

of r. The three unknowns T.J., M, and p are then obtained in the 
same manner outlined for the analysis method. ) 


Final Equations 

Auxiliary differential equation . - An auxiliary differential equa- 
tion for the pressure p in terms of T.j., M, and p is obtained from 
the equilibrium equations, which, after ejqiressing q^ and q^ in 
terms of q and p, combine to give 
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g dp 


1 dq^ C£ sec 3 


pq 


dr 2q^ dr h sin a 


(13a) 


But, 

£S! = 2:23! . ,pmZ 

g rgR*T 

so that, from equations (6b), (l3a) , and (l3b) 


(13b) 


1 dP 
P dP 




1 dM^ 


a + ^ M^/ dR 


1 dl4 


Tt 


2^ 


H cos 3 


(13c) 


where 


sin 


ct 


(13d) 


and where 


and 


P = 


R 


H = ^ 


hjp 


= H(B) 


~\ 


J 


(I3e) 


where is the compressor- inlet stagnation pressure^ is the 

impeller tip radius^ and h^ is the effective diffuser height at the 
impeller tip. Equation (l3c) is an auxiliary differential equation that 
relates the change in P to the change in and with radius E. 

Total temperature . - The change in T^ is given by equation (lie) 
or ( 12 ) ^ whi ch from e quat i ons (l3d) and (l3e) become 
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SR pq^HCp sin a 



and, for the case of Eeynolds* analogy. 


Tt ® 


t 

H cos P 



(14a) 


(14b) 


Mach number . - In order to determine the differential equation for 
the Mach number squared it is first necessary to express the second term 
of the continuity equation (?) in terms of known variables. From the 
meridional equilibrium equation (8) together with equations (4) and 
(13b) 

2 2 

1 dq^ tan P sec p 1 dp c^ sec p 

— _ - - _ (l5a) 

q^ dr r yM p dr h sin a 


The first term of the continuity equation (?) is expressed in terms of 
known variables by the equation of state (lOb) together with equation (6a 


1 dp 1 dp 

p dr p dr 



Tt dr 



1 dM^ 
M^ dr 


(15b) 


Substituting equations (l5a) and (l5b) into the continuity equation (?) 
and combining it with equations (l3d) and (l3e) result in 


1 dP 
P dR 



yM^ - sec^p 



1 dM^ 
M^ dE 


Tt dE 


^ 1 dH sec^ p 

+ + 

E 


H cos p 


H dE 
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1 dp 

which, combined with equation (l3c) to eliminate — — ^ finally gives 

P dR 


1 

dB - 


2o 

sec p 


? 2 1 

(l + yW - tan p) + 


2T^ dr 


{yl? - tan^p) ^ 


1 dE sec p 


H cos p H dE 


E 


(15c) 


Equation (l5c) determines the change in vith radius E along the 


1 dT4 


mean surface of revolution in terms of — 
equations (l4a) or (l4b). 


Tt dE 


which is known from 


Flow direction . - The differential of the flow direction p is 
obtained from equation (4a) 


1 d tan p 


tan p dr 


— ^ 


m 


dr 


which from the tangential equilibrium equation (9) and equation (l5a) 
becomes 


2 2 

1 d tan p sec p 1 dP sec p 

tan P dE P dE E 

and from equations (l3c) and (l5c) 



Equations (l4a) or (l4b), (l5c), and (16) are three differential equa- 
tions that can be solved simultaneously for and p. 
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Pressure . - After the variations in and 3 with radius 
E are known the pressure P can he obtained from the continuity equa- 
tion as follows? 


cos 3^ r^^ = pq cos 3 rh 


(17a) 


where the subscript 1 refers to known conditions (appendix C) at the 
diffuser inlet (where E is equal to 1.0 ) . Prom the equation of state 
(lOa) and from the definition of Mach number (equation (5)), equa- 
tion (l7a) becomes 


- COS 3^ 




Finally, from equations (6) and (l3e) 



1 cos 3^ 

EH cos 3 M 



(17b) 


Equation (l7b) determines P from the known conditions at the diffuser 
inlet and from the known values of T. , and 3 determined by the 

simultaneous solution of equations (ita) or (l4b), (l5c), and (l6). 

After the quantities P, '^t^ ^ known, all other 

quantities (p, T, q, qjj^, and q^) can be determined directly from 
equations (4), (s), (6), and the equation of state (lOa). 

Flow path . - The flow path on the mean surface of revolution in the 
vaneless diffuser can be obtained from the known variation in tan 3 
with E given by the solution of equation (16). From figure 7, which 
shows the flow path on a developed view of the mean siirface in the 
vicinity of E, 


tan 3 


B d0 
dR/sln a 


or 


d0 _ tan 3 
dE E sin a 


(18) 


Because 3 and a are known functions of E, equation (18) determines 
the flow path on the mean surface of revolution. 
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Influence coefficients . - In some analysis problems it may be con- 
venient or desirable to solve directly for one or more of the other 
dependent quantities rather than y? , and p. Also, in the design 

problem it may be desired to specify one of these quantities as a func- 

1 dH 

tion of E and to solve for the required value of — — . For these 

H dR 




T, q 

along the mean surface of revolution, as v^l as the 


IS 


and 


^0 


cases the change in the dependent variables 
with radius E 

change in T^, M^, and 3, must be expressed in terms of the known 

1 dT-i- t 1 

quantities ^ (given by equations (l4a) or (l4b)), g '' gog p -> 
t 


and — — (not known in the design problem), which quantities are 
H dR 


multiplied by influence coefficients. Thus, if X is any one of the 
dependent variables. 




dT^ 


I t T i T 1 

2 H cos 3 3 H dE 4 E 


(19) 


where through are influence coefficients that are determined 

in the same way that equations {l5c) and (l6) were developed. The 
influence coefficients for various dependent variables X are given in 
the following table: 


Influence Coefficients 


X 

I 

1 

I 

2 

^3 

^4 

p 

2 

7M 

2 1 2l 

7M [1+(7-1)MJ 

2 

-7M 

2 2 
-7M sec 3 

p 

sec 3 

M^(7 sec^3-tan^3) 

-M^ 

2- 

-M sec 3 

T 

7M^- sec^3 

(7-l)M^(7M^-tan^3) 

-(7-1)M^ 

- (7-1)M^ sec^3 


tan^3-l-7M^ 

2 (l-f^M^) (tan^3 - 7 ^) 

2(1+I^M^) 

2(1+ ^ M^)sec^3 

2 

q 

-2 

2(tan^P -yM^) 

2 

2 sec^3 

q 

m 

-sec^3 

M^(tan^3 - 7 sec^3) 

sec^3 

2„ .,2, 2_ 
sec 3 +M tan 3 

q0 

0 

sec^P - 

0 

sec^3 - 

tan p 

2„ 
sec 3 

sec^3 1+(7“1)M^ 

-sec^3 

2 2 

-M sec 3 
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If is equal to sec^JB^ the left side of equation (l9) is 

equal to zero because i ^ cannot be infinite. For this condition 

,2 2 ^ ® 

(M = sec p), which occurs when q^^^ is equal to the local speed of 
sound, that is when choke flow occurs (appendix D), equation (l9) 
becomes 




(l + sec^p) — + [i + (7-I) 

2 T. dR L 


sec p 


t 


H cos P 


1 

H dR 


sec^P 


= 0 
(19a) 


where the Influence coefficients given in the table have been substi- 
tuted in equation (19) with equal to sec^p. (All sets of influ- 

ence coefficients result in the same equation (19a).) Equation (l9a) is 
a condition that must be satisfied at the location of choke flow in a 
vaneless diffuser. In particular, if heat transfer and friction are 
absent, equation (l9a) becomes 


_1_ d<BE) -tan^p 

HR dR E 


Because HR is directly proportional to the flow area, it is seen that 
choke flow does not occur at the throat, or position of minimum flow area 

^ ^ ^ = 0 jj, but at a point where the flow area is decreasing in the 
increasing E. 


dR 

direction of 


Small-stage efficiency . - The small-stage, or polytropic, efficiency 
at a given radius E on the mean surface of revolution in a vaneless 
diffuser is defined as the ratio of the ideal (ignoring friction and heat 
transfer) to the actual differential change in static enthalpy with 
radius required to accomplish the actual differential change in static 
pressure with radius. This definition leads to the following expression 
for the small-stage efficiency (appendix E) 


n = 


1 ^ 
P dR 


1 ^ 
P dR 





-i- ^ 

dE H cos P 


(20a) 


Equation (20a) indicates that in the absence of heat transfer 
(dT^/dR = 0) and friction (^ = 0) the small-stage efficiency is 100 
percent. Also, for heat transfer from the fluid to the diffuser walls. 
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dT. 


1 

— — — is negative and therefore results in an apparent increase in the 


Tt dE 


small- stage efficiency, 
small- stage efficiency. 


Thus, in the presence of heat transfer, the 
as just defined, is not a good measure of the 
performance of vaneless diffusers in that it is not a measure of the^^ 
magnitude of the losses involved. In the absence of heat transfer 

/ / X ® 

is zero and equation (20a) can be rearranged to give (appendix E) 


T1 = 1 


t (M^ - sec^p) 

t (rM^ - tan^p) - cos 3 ^ 


(20b) 


If is equal to sec^p (choke flow condition, appendix D), the 

efficiency given by equation (20b) becomes indeterminate, because the 
numerator of the fraction is zero and, from equation (l9a), so is the 
denominator . 


NUMERICAL PEOCEDUEE 

A specific numerical procedure is outlined for both the analysis 
and design problems, however, any other standard numerical procedure can 
be used. In the analysis problem the variation in fluid properties 
with E are determined for a specified geometry of the vaneless dif- 
fuser. In the design problem the variation with E in one of the fluid 
properties is prescribed and the remaining fluid properties together 
with the variation in diffuser height H with radius E are deter- 
mined. The numerical procedures for both the analysis and design prob- 
lems are essentially the same but differ in detail and are therefore 
discussed separately in this section. 


Analysis Problem 

Primary quantities . - In the analvsis problem the variation with 
E in three primary quantities (T^, and p) are obtained from 

three differential equations; (l4a) or (l4b), (l5c), and (l6). These 
equations are nonlinear and it is necessary to solve them by numerical 
methods. The suggested stepwise procedure is as follows: 

(l) The values of T^, M^, and p at the diffuser inlet are 

estimated from the impeller design and operating characteristics 
(appendix C). 
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(2) At any radius E, if y? , and tan 3 are known, the 

change in (that is, AT-j^) for a small increase in E (^that is, 

AE) is computed directly from equation (l4a) or (l4b). 

(3) The change in (that is, AM^) for the same small change 

in E (that is, AE) is computed directly from equation (l5c) in which 

1 ®t / \ 

- — — — is obtained from step (2). 

T^ dR 

(4) Finally Atan P is computed from equation (16). 


(s) Thus, at (e + AE) the approximate values of T-^, , and 

tan p are known from the values of T-t, M^, and tan p at E and 
the approximate values of AT^, AM^, and Atan 3 given by steps (2) 
through (4). 


(6) At 

1 d tan P 
tan p dE 
and (l6) from 
in step (5 ) . 


(e + AE) approximate values of ~1_, and 

Tt dR dR 

can then be determined by equations (l4a) or (l4b), (l5c) 
the approximate values of T-j^, M^, and tan p obtained 


(7) The 

(E + AE) are 

and so forth, 
and (6), 


final values of AT-^, AM^, and Atan p between E 
obtained directly from the arithmetic average of i 
at E and (e + AE) as given in steps (2), (3), (4^', 


and 
dTt 
dR ^ 


(8) The values of T^, 1^, and tan p at (E + AE) are deter- 
mined from the known values of T-j^, M^, and tan P at E and from 
the values of AT^, AM^, and Atan P obtained in step (7). 


(9) The stepwise procedure outlined in steps (2) through (8) is 
repeated for small values of AE starting at E equal to 1.0 (where 
T. , M^, and tan 3 are obtained by step (1)) and continuing to the 

diffuser exit. For the numerical examples of this report AE was 0.02, 
0.03, and 0.05 for the first three increments and 0.10 for the remaining 


Secondary quantities . - After the distribution of T^, M^, and 

tan p with E are known, the distribution of P, p, T, q, qjjj, and q^ 
can be determined directly from equation (l7b) and from equation (4), 

( 5 ), (6), and the equation of state (lOa) . 


Flow path . - The flow path on the mean surface of revolution in the 
vaneless diffuser is given by 0 as a function of E along the sur- 
face. Because tan 3 and sin a are known functions of E, the flow 
path (9 = 9 (E)) can be determined by the numerical integration of 
equation (I8) assuming 9 equals zero at E equals 1.0. 
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Efficiency. - The small-stage^ or polytropic, efficiency t] at 
each radius E is given by equation (20a). The diffuser efficiency 

between radius B equals 1.0 and B is given by equation (E6) 
in appendix E. 

Design Problem 

In the design method the variation in effective diffuser wall 
spacing with radius is determined for a prescribed variation in one 
fluid property. Eor efficient diffuser designs the selection of the one 
fluid property and its optimum prescribed variation will depend on 
viscous flow effects that are considered in boundary- layer studies but 
will not be investigated in this report. (Nor is the magnitude of the 
boundary- layer displacement, required to obtain the geometric wall 
spacing from the effective spacing, considered in this report.) 

Auxiliary equation . - In the design problem the variation in H 
with E is unknown and must be determined to satisfy a specified varia- 
tion in one characteristic of the flow (qjjj, for example) with E. From 

1 dH 

this specified variation in one characteristic of the flow — — can he 

determined from equation (19). The quantity AH between E and 

1 dH 

(E + AE) is obtained from the average value of -5^ at E and 

JbL CLn 

(E + AE) . After H has been determined at (E + AE) the final values 
of AT.^^ AM^^ AP^ and Atan p between E and (E + Z^) are 
obtained by the same procedure previously outlined for the analysis 
problem. The process starts at E equals 1.0 (where H equals 1.0) 
and is repeated for specified increments of E up to the diffuser exit. 

Com.plete solution . - After the variation in H^ T^^ and 

tan P with E are known, the variation in P, T, q, q^, and q^ 

can be determined directly from equation (l7b) and from equations (4;, 
(5), (6), and the equation of state (lOa). The flow path is determined 
by equation (l8) as outlined previously, and the small-stage, or poly- 
tropic, efficiency is determined by equation (20a). 


NUMERICAL EXAMPLES 

The numerical examples of this report are divided into three groups: 
(l) effects of some operating conditions, (2) effects of diffuser wall 
spacing, and (3) a vaneless diffuser design problem. 
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Effects of Some Operating Conditions 

The first group of numerical examples shows the effects of heat 
transfer and friction on the flow in vaneless diffusers. Three numer- 
ical examples are given: {l) isentropic compressible flow^ (2) com- 

pressible flow with friction^ and (3) compressible flow with friction 
and heat transfer. 


Inlet conditions . - Eor the first group of numerical examples the 
flow conditions at the diffuser inlet (E = l.O) are: 


= 3.022 
= 1.370 
(Tt)i = 941° E 


(tan = 3.829 

These conditions were estimated (by methods given in appendix C) for 
the following design and operating conditions of the impeller: 


Compressor flow coefficient, cp 0.75 

Impeller tip Mach number, Mji 1.5 

Impeller slip factor, p 0.9 

Impeller polytropic efficiency, t] 0.9 

Compressor stagnation inlet temperature, T^, °E 520 


Diffuser design . - Thfe design characteristics of the diffuser are: 


Passage height, H 

(constant flow area normal to q^j^) 

Wall temperature, T^, °R 

Friction parameter, t 

Cp = 0.003 (a relatively low value, see reference 4) 
=10 (so that sin a is constant) 


1 / ^T 

sin 



750 

0.030 


The Eeynolds* analogy was used to determine the heat-transfer coefficient 
so that equation (l2) was used to determine the change in total tempera- 
ture with radius. 


Results . - The results of the first group of three numerical 
examples are given in figure 8. In figure 8(a) is shown the change in 
W- with E for the three numerical examples. The effect of friction 
is to reduce at each radius E (because, although the smaller 
meridional velocity component q^ is increased as usual, the larger 
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tangential component is decreased) and the effect of heat transfer 

from the fluid is to increase slightly (primarily because of the 

reduced speed of sound at the lower temperature) for the magnitudes of 
and T-y involved in these examples. 

In figure 8(b) is shown the change in P with R. As expected, 
the effect of friction is to reduce P at each radius (primarily 
because of the decreased values of which require a smaller pressirre 

gradient for equilibrium) . The effect of heat transfer from the fluid 
is to raise P slightly for the magnitudes of T-j. and involved in 

these examples. 

In figure 8(c) is shown the change in flow direction 6 with E. 
The effect of friction is to reduce (3 because qg is reduced and q^ 
is increased to satisfy continuity with lower density due to lower P. 
The effect of heat transfer from the fluid is to increase p slightly 
because of the reduced value of q resulting from the increased value 
of p. 

In figure 8(d) is shown the flow path in the vaneless diffuser. 

The effect of friction is to shorten the flow path because p is 
decreased (fig. 8(c)). The effect of heat transfer is to lengthen the 
path slightly. 


In figure 8(e) is shown the change in polytropic, or small-stage, 

efficiency t) with radius E. The effect of friction is to reduce the 

efficiency at each radius. The effect of heat transfer from the fluid 

is to increase t) greatly. In fact for the larger values of E where 

2 1 ®t 

M is relatively small the term involving — becomes greater than 

t 

the term involving the friction parameter t (see equation (20a)) and. 


because for heat transfer from the fluid — 


dT+ 


T+ dR 


is negative, tj 


becomes greater than 100 percent. 


For the example with friction but no heat transfer it is interesting 
to note that, although the friction losses must be greater at the lower 
values of E because of the larger velocities, the polytropic effi- 
ciency is higher. From equation (20a) the higher efficiency must result 

1 dP 

from a higher rate of pressure rise "™ ^ i+v 


— — compared with 
P qP 


H cos p 


at 


the lower values of E. In current compressor designs the local poly- 
tropic efficiency at the lower values of E will be considerably reduced 
(reference 4, for example) because of mixing losses resulting from the 
nonuniform flow conditions at the impeller discharge. 
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The diffuser efficiency .(^D^2 ^ equal to 2.0 is determined hy 

equation (E6) of appendix E and is indicated for each of the three numer- 
ical examples in figure 8(e). The value of ('n 7 )) 2 ^ "the example 

with heat transfer, is considerably less than the values of the poly- 
tropic efficiency t] at the larger values of E for this example, 
because for these larger values of E the rate of pressure rise is 
smaller and therefore tj has less effect upon the value of ^Ap) 2 * 

Even so, the value of ('np )2 example with friction and heat 

transfer is almost 100 percent. But from figure 8(b) the pi^ssure P 
at E equal to 2.0 is not much different for the two examples with and 
without heat transfer (but with friction) so that the losses are about 
equal for these examples and (ilj)) 2 ^ which is considerably different, 
is therefore not a reliable measure of the losses when heat transfer 
effects are present. 

For the example with friction and no heat transfer the value of 
(^1)^2 0.824 in spite of the relatively low friction 

coefficient (c^ = 0.003) and in spite of neglecting the mixing losses 
due to nonuniform flow conditions at the impeller discharge. Thus, the 
friction losses in most vaneless diffuser designs are considerable and 
result from the (usually) large ratio of wetted surface to flow area. 

The diffuser efficiency can be improved by lower values for S (pro- 
vided other design and flow variables remain unchanged) and these lower 

values for ^ can result from lower values of v— , which means, for 

hp-’ ’ 

example, larger compressor flow rates for a given impeller tip radius. 

A general conclusion resulting from the first group of numerical 
examples is that heat transfer from the fluid has the opposite effect of 
friction on pressure rise in vaneless diffusers and is therefore to be 
desired. Heat transfer to the fluid, on the other hand, can be expected 
to have the same effect as friction and is therefore to be avoided. 


Effects of Diffuser Wall Spacing 

The second group of numerical examples shows the effects of passage 
height h (that is, spacing of the diffuser walls normal to the mean 
surface of revolution) on the flow in vaneless diffusers. The losses 
in a vaneless diffuser should increase with the velocity squared, with 
the ratio of wetted perimeter (at each radius) to diffuser wall spacing 
(that is, with the ratio of friction area to flow area), and with the 
length of the flow path in the vaneless diffuser. For a given com- 
pressor flow rate the square of the velocity and the ratio of wetted 
perimeter to diffuser wall spacing increases as the diffuser wall spacing 
h is decreased, but the length of the flow path decreases. The object 
of the second group of numerical examples is to determine the relative 
magnitudes of these opposing effects on the losses in vaneless diffusers, 
and to determine the optimum wall spacing h, if such an optimum exists. 
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Diffuser design . - As for the first group of numerical examples the 
diffusers of the second group were designed for constant flow area 

(H = E“^) but with the diffuser wall spacing hj at the impeller dis- 
charge (diffuser inlet), and therefore throughout the diffuser, varying 
(among examples) over a wide range. The values of bj were selected 
for the eight numerical examples of the second group such that the 
friction parameter ^ varied from 0.010 to 0.038 in seven increments of 
0.004. (^ varies inversely with hj, equation (l3d).) The friction 

coefficient was assumed to be constant, that is independent of diffuser 
wall spacing, and therefore the possibility of separated flow for large 
spacing of the diffuser walls was not considered. Heat-transfer effects 
were not considered. The example for t, = 0.030 was the same as that in 
the first group of numerical examples with friction but no heat transfer. 

Inlet conditions . - For the second group of numerical examples the 
diffuser inlet conditions varied with the diffuser wall spacing because 
for a constant compressor flow rate W the compressor flow coefficient 
cp (equation (Cll), appendix C) varies inversely with the passage height 
hj at the impeller tip. (Note that the blade height at the impeller 
tip is also assumed to vary with hj. and this variation is assumed to 
have no effect on the impeller efficiency, and so forth. ) For the 
selected variation in hp the flow coefficient varies from 0.25 to 0.95 
in seven increments of 0.10. The remaining design and operating condi- 
tions of the impeller are the same as for the first group of examples. 
Thus, 


Impeller tip Mach number, Mp 1.5 

Impeller slip factor, |.i 0.9 

Impeller polytropic efficiency, 0.9 

Compressor stagnation inlet temperature, T^, 520 


The flow conditions at the diffuser inlet were estimated (by methods 
given in appendix C) from the impeller design and operating conditions 
and are given in the following table; 


Example 


cp 

Pi 

Mi^ 

(tan p)^ 

(°R) 

a 

0.010 

0.25 

3.174 

1.272 

11.879 

941 

b 

.014 

.35 

3.157 

1.283 

8.453 

941 

c 

.018 

.45 

3.133 

1.298 

6.541 

941 

d 

.022 

.55 

3.103 

1.317 

5.317 

941 

e 

.026 

. 65 

3.066 

1.341 

4.462 

941 

f 

.030 

.75 

3.022 

1.370 

3.829 

941 

g 

.034 

.85 

2.970 

1.406 

3.339 

941 

h 

.038 

.95 

2.909 

1.448 

2.945 

941 


d 
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Eesults . “ The results of the second group of eight numerical 
examples are given in figure 9. All results are based on the assumption 
that as the diffuser wall spacing increases^ that is as cp descreases^ 
the friction coefficient remains constant and flow separation does not 
occur on the diffuser walls. In figure 9(a) is shown the change in tT 
with R for the eight examples. The effect of decreasing cp (that is, 
of increasing hj) is to decrease at each value of E. This 

decrease in results primarily from the decrease in qj^ resulting 

from the increased flow area that occurs when hp is increased. 

In figure 9(b) is shown the change in P with E for various 
values of cp. The effect of decreasing cp is to increase P because 
is decreased (fig. 9(a)). 

The change in p with E for various values of cp is shown in 
figure 9(c). As cp is decreased the velocity cotaponent q decreases 
so that the flow direction p increases as shown. As cp approaches 
zero, p approaches 90° for all values of E. 

In figure 9(d) is shown the flow path in the vaneless diffuser for 
the various values of cp. The effect of Increasing cp is to decrease 
the length of the flow path because p is reduced (fig. 9(c)). 

The change in polytropic efficiency t) with E for the various 
values of cp is shown in figure 9(e). The effect of decreasing cp is 
to increase ti at the larger values of E. As the value of cp 
decreases the length of the flow path increases (fig. 9(d)), which 
should increase the diffuser losses, but (fig. 9(a)) and the fric- 

tion parameter t decrease, which should decrease the losses. At cp 
equal to zero P is 90°, or cos p is zero, and ^ is zero so that 
the ratio l/cos p contained in the expression for q (equation (20a)) 
becomes indeterminate. However, extrapolation of the results in 
fis^r® 9(f) indicates that q has its peak value for <p equal to zeroj 
and thus, if separation does not occur and if the friction coefficient 
Cj is unaffected by the diffuser wall spacing, the diffuser efficiency 
is always Improved (slightly, see fig. 9(f)) by spacing the diffuser 
walls farther apart. Furthermore, if the diffuser walls are spaced 
farther apart and the compressor flow rate is increased proportionately 
so that cp remains unchanged, the diffuser efficiency should be Improved 
markedly because the ratio of wetted surface to flow area is decreased 
without Increasing the length of the flow path in the diffuser and with- 
out introducing the risk of boundary- layer separation, which must other- 
wise be expected if the diffuser walls are spaced far apart. Thus, as 
also concluded from the first group of numerical examples, the diffuser 
efficiency can be improved by Increased compressor flow rates for a 
given impeller tip radius so that the diffuser walls can be spaced 
farther apart without resulting in boundary- layer separation or 
Increasing the length of the flow path. 
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In figure 9(f) is shown the effect of cp on the diffuser efficiency 
As expected from figure 9(e), (’lj ))2 increases as* cp decreases, 
but the rate of increase is less for the smaller values of cp and in no 
case are large gains in efficiency to be realized by decreasing cp, 
that is, increasing the diffuser wall spacing. Thus, unless the flow 
rate is increased proportionately so that cp remains constant, very wide 
spacing of the diffuser walls is not recommended because of (only) a 
small gain in efficiency and a great risk of boundary- layer separation. 


The third part of the section on numerical examples is a sample 
vaneless diffuser design problem. The design variables in a vaneless 
diffuser are 


equal to 90°) and the design problem will be to determine H(E) for a 
prescribed variation in q^. 

For purposes of demonstrating the design method it is assumed that 
the deceleration of q^ is the. criterion for boundary- layer separation 
in a vaneless diffuser and that the criterion is that given in refer- 
ence 5, page 159, so that a safe rate of deceleration is 


where 6 is proportional to the boundary- layer thickness. For purposes 
of this design example it is assumed that 6 is equal to h/2, which is 
the effective thickness of a fully developed boundary layer in the vane- 
less diffuser. Thus, 


A Vaneless Diffuser Design Problem 


H = H(E) 




or 
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1 ^ ^ 
q dR ~ H 

if rip/h^ is equal to 10. Because of the assumptions involved this 

specified variation in with H may have no practical significance 

with regard to vaneless diffuser performance and has been selected only 
to demonstrate an application of the design method. It should be 
pointed out that design variations In H affect primarily the velocity 
component qj^ and through this component the flow direction 3. 

Inlet conditions . - The impeller design and operating conditions 
are the same as for the first group of numerical examples and so the 
diffuser inlet conditions are the same 


= 3.022 
= 1.370 

(T^)^ = 941° E 
(tan 3)^ = 3.829 

Diffuser design . - The variation in H with E is to be deter- 
mined. Heat-transfer effects are neglected, and the value of the fric- 
tion parameter t is the same as for the first group of numerical 
examples (0.030). 


Ee suits . - The results of the design problem are given in figure 10. 


In figure 10(a) is shown the variations in H, — 


T) With radius E. As specified, — 

' <im ® 


<lm 


is equal to -H 


r-1 


In order 


and 


to accomplish this variation, H at first decreases with increasing E 
and then increases to approximately its initial value at E equal to 
2.0. At the larger values of E this variation in H results in some- 
what wider spacing of the diffuser walls than existed in the previous 
numerical examples where the wall spacing decreased continuously with 
increasing E in order to maintain a constant flow area normal to q 
As a result of this wider spacing of the diffuser walls the values of^’u 
are somewhat higher (in keeping with the results of the second group of 
numerical examples) than for the previous examples with the same values 
of cp and i. The variation in 3 with E was slightly more than 3° 
so that the flow path (fig. 10(b)) is approximately a logarithmic spiral. 
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SUMMAEI OF EESULTS AND CONCLUSIONS 

An analysis method and a design method have been developed for one- 
dimensional, compressible flow with friction, heat transfer, and arbi- 
trary area change in vaneless diffusers with arbitrary profiles in the 
axial-radial plane. The effects of mixing losses due to nonuniform flow 
conditions at the impeller discharge are not considered. In the analysis 
method the variation in fluid properties, including the velocity and 
flow direction, can be determined as a function of radius for a prescribed 
variation in diffuser wall spacing with radius. In the design method 
the variation in diffuser wall spacing and all fluid properties except 
one can be determined as a function of radius for a prescribed variation 
in the one fluid property. For efficient diffuser designs the selection 
of the one fluid property and its optimum prescribed variation will 
depend on viscous flow effects that are considered in boundary- layer 
studies but are not investigated in this report. 

Three groups of numerical examples are presented in which the 
effects of friction, heat transfer, and diffuser wall spacing are 
investigated; and a sample design problem is presented. As a result of 
these examples it is concluded that: 

1. Heat transfer from the fluid has the opposite effect of friction 
on pressure rise in vaneless diffusers and is therefore to be desired. 
Conversely, heat transfer to the fluid has the same effect as friction 
and is therefore to be avoided. 

2. If the friction coefficient is unaffected by the diffuser wall 
spacing, and if flow separation does not occur, the diffuser efficiency 
is Improved slightly (for a given compressor flow rate) by spacing the 
diffuser walls farther apart. 

3. Even with relatively low friction coefficients and neglecting 
mixing losses at the impeller tip, the friction losses in most vaneless 
diffuser designs are considerable, as indicated by computed diffuser 
efficiencies in the low 80* s, and these losses result from the usually 
large ratio of wetted surface to flow area in vaneless diffusers. 

Diffuser efficiencies can be improved by increased compressor 
flow rates for a given Impeller tip radius so that the diffuser walls can 
be spaced farther apart (thus, reducing the ratio of wetted surface to 
flow area) without increasing the length of the flow path in the diffuser. 

5 . In the presence of even small heat- transfer effects the usual 
definition of diffuser efficiency, whleh definition neglects corrections 
for heat transfer, is not a measure of the diffuser losses. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 11, 1951 
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APPENDIX A 
SYMBOLS 


3^ 


Cf 


Cp 

Ct 

g 

H 


h 

h* 


i 

J 

M 

Mrr 


n 

P 

P 

Po 

dQ 


The following symbols are used in this report: 

annulus flow area at impeller tip^ equation (C12) 
local speed of sound 

skin-friction coefficient^ equation (B2) 
stagnation speed of sound upstream of impeller 
specific heat at constant pressure 
local stagnation speed of sound 
acceleration due to gravity 

effective passage height ^ or diffuser wall spacing, 
ratio, h/hrp 

effective passage height, or diffuser wall spacing 
coefficient of heat transfer, equation (lla) 

I^, I^, influence coefficients, equation (l9) 

enthalpy 

mechanical equivalent of heat 
local Mach number, equation (S) 
impeller tip Mach number, equation (C3) 
polytropic exponent, equation (C6) 
pressure ratio, p/Pq 
static (stream) pressure 
stagnation pressure upstream of impeller 
heat transfer rate from the fluid 
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velocity 

'3m 

velocity component in meridional^ or axial-radial^ 
plane ^ eq-uation (3) 

*3-r ^ ^0 •> 'Iz 

velocity components in r-^ 6-^ and z-directions^ 
respectively 

R 

radius ratio ^ r/r^ 

R* 

perfect gas constant 

r, e, z 

cylindrical coordinates^ 0 considered positive in 
counterclockwise direction when r0-plane is viewed 
from negative z-direction 

T 

static (stream) temperature 

To 

stagnation temperature upstream of impeller 

Tt 

local stagnation^ or totals temperature 

T 

■^w 

diffuser wall temperature 

t 

time 

W 

compressor flow rate 

X 

dependent variable 

a 

slope of center line between diffuser walls in meri- 
dional^ or axial-radial plane ^ equation (2) 

P 

flow direction on mean surface of revolution between 
diffuser walls ^ equation (4a) 

r 

ratio of specific heats 

A 

small finite increment 


friction parameter^ equation (l3d) 


polytropic ^ or small- stage ^ efficiency^ equation (E1) 

(■^D^R 

diffuser adiabatic efficiency based upon change in flow 
conditions for change in radii from 1.0 to R, equa- 
tion (E5) 
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4 

P 

Po 

Pt 

T 

cp 

O) 

Subscripts ; 

a 

i 

R 

T 

1 , 2 


impeller slip factor 

static (stream) wei^t density- 

stagnation density upstream of Impeller 

local stagnation density 

shear stress due to skin friction 

compressor flow coefficient, equation (Cll) 

angilLar -velocity of lng)eller 

actual 

ideal 

value at R 
impeller tip 

value at R equal to 1,0 or 2.0 
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APEEKDIX B 


EQUILIBRIUM EQUATIONS 

Meridional and tangential equilibrium equations are developed for a 

fluid particle on the mean surface of revolution in a vaneless diffuser. 

% 

Meridional equilibrium. - The equation for meridional equilibrium 
of a fluid particle ^fig. 6) in the direction of on the mean surface 

of revolution is obtained from a balance of the pressure and shear 
forces with the force required for acceleration. 

The differential pressure forces (opposed to the direction of 
are equal to the differential change of end forces on the particle minus 
the component of the differential side forces on the particle in the 
direction of q^j^, 

„ . . . ^ d(phr d0)dr pd(hr d0)dr 

Differential pressure forces = ^ (Bl) 


where the component of the differential side forces in the direction of 
qjj^ (last term in equation (Bl)) is equal to the pressure p multiplied 
by the projected area (in the direction of q^) of the side surfaces of 
the particle (fig. 6) . 


The differential shear stress T on a diffuser wall is 
the direction of q and is given by 



opposed to 


(B2) 


where c^ is the skin friction coefficient. The differential shear 
forces in the meridional direction on the fluid particle in figure 6 are 
opposed to the direction of qj^ and act on both walls of the diffuser. 
From equation (B2) , 

^ r d0 dr 

Differential shear forces = 2 t cos p — : 




pqi 


cos P 


r d0 dr 
sin a 


(B3) 


The acceleration of the fluid particle in the direction opposed to 
qjQ is made up of; (l) the component of the centripetal acceleration 


iqjji 

sin a, and (2) the negative of the acceleration But, 
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^ dr 

dt “ dr dt “ ^ dr 


so that the differential force required for acceleration of the fluid 
particle in figure 6 becomes 


/Differential force requiredX 
for acceleration in direc- I = ^ : 


ytion opposed to 




g sin a 



- 


d5^ sin 


(B4) 

The sum of the differential pressure forces and shear forces must 
equal the force required for acceleration so that from equations (Bl) , 
(B3), and (B4) 

g dp , P 1 q^ ^ 

p dr h sin a r “ dr 


Equation (8) is the equation for meridional equilibrium of a fluid par- 
ticle on the mean surface of revolution in a vaneless diffuser. 

Tangential equilibrium. - The equation for equilibrium of a fluid 
particle (fig. 6) in the tangential direction on the mean surface of 
^ revolution is obtained from a balance of the shear forces with the force 
required for acceleration. 


The differential shear forces in the tangential direction on the 
fluid particle in figure 6 are opposed to the direction of q^ and act 
on both walls of the diffuser. From equation (B2) ^ ” 

Differential shear forces = 2t sin 6 ^ 

sin a 


„ pq . „ r d0 dr 
= c^ sin p ^ 


g 


sin a 


(B5) 


The tangential acceleration of the fluid particle opposed to the 
direction of q^ is made up of: (l) the negative of the tangential 


acceleration r ^•nd (2) the negative of the Coriolis accelera- 

. . „ s in a 

tion 2qj^qg 


■e 

r y 
But, 
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dt 


% dr 


r dt 


,2 dt 


1 *^<16 dr 'le dr 


r dr dt 


,2 dt 


r \dr r 


so that the differential force required ior acceleration of the fluid 
particle in figure 6 becomes 


( Differential force requiredX , (S 

for acceleration in direc- | ts - £ — ■ 
tion opposed to q 


■e 




sin a 


<100 <3^<le sin a 

q^^sina + ^ 


(B6) 


The differential shear force must equal the differential force 
required for acceleration so that from equations (B5) and (B6) 


c^q^ sin p 


'^'^6 •Vh 

h sin a dr r 


( 9 ) 


Equation (9) is the equation for tangential equilibrium of a fluid 
particle (fig. 6) in the tangential direction on the mean surface of 
revolution in a vaneless diffuser. 
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APPEEDIX C 

ESTIMAIED VALUES OF P^ AND TAN P AT DIFFUSER 

INLET (IMPELLER TIP) 

Total temperature. - The total ten5)erature (T^)]_ at the diffuser 

inlet, or impeller tip, can be obtained from the steady flow energy equa- 
tion, where for convenience heat-transfer effects have been considered 
negligible, 

+ 2^ (Cl) 

where p is the impeller slip factor and cjo is the angular velocity of 
the impeller so that p(aiTm)^/g is the impeller work per pound of fluid. 
But 

P v-1 

and 


(%)l = 


so that dividing equation (Cl) by JCpT^ 


— = 1 + Ili 

To ^+2 


(2p - p^) 


('3m); 




vhere the impeller tip Mach number jy^ is defined by 


CDT, 


and 


(le) 


~ (equal to p^Hj,^) are equal to zero so that 


(C2) 


Kp = — 

The total temperature (T^)3_ is given by equation (C2) when 




(C3) 

2 


(^t)i = To[i + (r- 1 ) pMj.2] 


(C4) 


Pressure. - The pressure p^ at the diffuser inlet is obtained 
from the temperature ratio (equation (C2)) by 
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!i-P -f!i 

Po ^ 


n 

n-1 


(C5) 


where the flow in the impeller^ which involves viscous losses, is rep- 
resented hy a polytropic process for which the polytropic exponent n 
is related to the polytropic efficiency of the impeller ri hy (refer- 
ence 8, p. 449, for example) 


n r 

n-1 " ^ y-1 


(C6) 


The quantity (q^)p in equation (C2) is unknown but will be determined 
later from continuity considerations. 

In like manner the density is related to the temperature ratio 

(equation (C2)) by 

1 




(C7) 


Mach number. - The local lyfech number squared (M )-| at the diffuser 
inlet is defined by 


(m2)^ = — ^ 


+ (q (q^)-,^ + (q T 


m 1 


e 1 


ml _o^ 

t7 


(C8) 


where Tq/T-]_ is given by equation (C2) , where 


Z “ ^T 


(C9) 


and where 
continuity 


(Vi 


is detei*mined from continuity considerations. From 


W = 2rtr^h^p3_ (qJi 


or 


1 _ cp 

“o ’ fl/Po 


(CIO) 
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where the compressor flow coefficient cp is defined by 


qp = 


W 


Po^^o 


(Cll) 


(C12) 


in which the annulus area at the impeller tip a^ is given by 

Equations (C7) and (CIO) are solved simultaneously for ^ so that 

equation (C8) can be solved for 


Flow direction. - -The tangent of the flow direction p-, at the 
diffuser inlet is defined by 


(tan p)^ = 


<V°o>: 

(V<=o): 


(C13) 


where and are given by equations (C9) and (CIO), 

respectively. 

2 

Thus^ Mq , and (tan p)j^ are estimated by equa- 

tions (C4), (C5), (C8), and (C13) , respectively. 
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APPENDIX D 


CONDITION FOR MAXIMUM, OR CHOKE, FLOW IN VANELESS DIFFUSERS 


If W is the flow rate through a vaneless diffuser 


W = 2nrhpq^ 


(Dl) 


where p is related to the stagnation density 
hy (reference 7, p. 26, for example) 


at a given radius r 


o> 

C\3 



1 

r-1 


(D2) 


(in the presence of heat transfer and friction p^ is a function of 
radius.) At any given radius as is increased from low values W 

increases until a maximum^ or choke^ flow occurs. This maximum occurs 
when 


or^ from equation (D1)^ 



(D3) 


but, from equation (D2) , 



(D4) 


where from reference 7, page 26, for example 



so that equation (D4) becomes 



(D5) 
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where 


After equations (D3) and (dS) are 

Oja = c (D6) 

so that the maximum^ or choke^ flow occurs in vaneless diffusers when 
the meridional component of velocity is equal to the local speed of 

sound c. Expressed in terms of equation (D6) becomes 

0 =s - sec^ p (D7) 

which is the condition for maximum, or choke, flow in vaneless diffusers. 



combined 
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APPENDIX E 


SMALL-STAGE EFFICIENCY AND DIFFUSER EFFICIENCY 


The small-stage, or polytropic, efficiency t] at a given radius 
R on the mean surface of revolution in a vaneless diffuser is defined 
as the ratio of the ideal to the actual differential change in static 
enthalpy with radius required to accomplish the actual differential 
change in static pressure with radius 


where the ideal differential change in static enthalpy with radius 
(di/dr) is given hy (reference 6, p. 102) 


and where the actual differential change in static enthalpy with radius 
( di/dr is hy definition 


Equation (El) is the usual definition of small-stage, or polytropic, 
efficiency and assimies that heat-transfer effects are negligible. From 
equations (El) to (E3) and equation (l3e) 



(El) 



(E2) 



(E3) 


1 dP 


Ti 


r 1 

r-1 T dR 


(E4) 


which from equations (6a) and (l3c) becomes 



1 ^ 
P dR 


(20a) 



+ 


H cos p 
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Equation (20a) gives the small-stage^ or polytropic, efficiency in terms 
of the local pressure differential, the total temperature differential 
and the parameter which Involves the local skin-friction coeffi- 

cient. 


1 d.T 

In the absence of heat transfer ^ equals zero and from equa- 

±t dK 


tlons (l3c) and (l5c) 


1 dP 

V - 




2 

sec p 


- tan^ p) 


_____ 1 ^ sec^ p 

H cos p ~ H dR " R 




] 


H cos p 


so that equation (20a) becoines 


T) = 1 - 


- sec^ p) 


^(tM^ - tan^ P) - cos P (^ + 


dH , H sec2 p 


(20b) 


Equation (20b) expresses ti in terms of the friction parameters ^ and 
the diffuser geometry. 


The diffuser efficiency (fl[))R^ which measures the diffuser per- 
formance between the diffuser inlet at R equals 1.0 and a point R 
on the mean surface of revolution in the vaneless diffuser, is defined 
as the ratio of the ideal 'to the actual static temperature rise required 
to accomplish the actual static pressure rise between the radii 1.0 and 
R. 




(Tr - Ti), 
- Tija 



(E5) 


where the ideal temperature ratio 



is related to the actual 



static pressure ratio by 
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so that equation (E5) becomes 




r-1 

r 



(E6) 


Equation (E6) gives (i1d)r in terms of the known values of P and T 
at the radii 1.0 and R. If R approaches 1.0, 


r-1 r-1 



r 


r 




-^1 


+ 


r-1 

r 



and 


so that 




(Vi -I 


1 ^ 

P dR 
r 1 dT 
vr-1 T dR/ 


which corresponds to the definition for the small-stage efficiency given 
by equation (E4) . 


REFERENCES 

.1. Polikovsky, V., and Nevelson, M. : The Performance of a Vaneless 

Diffuser Fan. NACA IM 1038, 1942. 

Brown, W. Byron, and Bradshaw, Guy R.: Method of Designing Vaneless 

Diffusers and Experimental Investigation of Certain Undetermined 
Parameters. NACA TN 1426, 1947. 


2 . 


2391 


2391 


6S 


MCA TN 2610 


41 


3. Shapiro, Ascher H. , and Hawthorne, W. R. : The Mechanics and Thermo- 

dynamics of Steady One-Dimensional Gas Flow. Jour. i^pl. Mech., 
vol. 14, no. 4, Dec. 1947, pp. A317-A336. 

4. Brown, W. Byron: Friction Coefficients in a Vaneless Diffuser. 

NACA TN 1311, 1947. 

5. Prandtl, L. : The Mechanics of Viscous Fluids. Vol. Ill of Aero- 

dynamic Theory, div. G, sec. 13, W. F. Din-and, ed., Julius Springer 
(Berlin), 1943. 

6. Keenan, Joseph H. : Thermodynamics. John Wiley & Sons, Inc. (New 

York), 1941. 

7. Liepmann, Hans Wolfgang, and Puckett, Allen E.; Introduction to 

Aerodynamics of a Compressible Fluid. John Wiley & Sons, Inc. 

(New York), 1947. 

8. Wislicenus, George F.: Fluid Mechanics of Turbomachinery. McGraw- 

Hill Book Co., Inc. (New York), 1947. 


42 


MCA TN 2610 



Figure 1. - Assembly of vaneless diffuser with mixed-flow Impeller. 
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Figure 3. - Diffuser profile, velocity components, and 
coordinates in meridional, or axial -radial, plane. 
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Fluid particle on surface of revolution 

diffuser walls. 


generated by center line between 


Figure 4. 
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Figure 5. -Mean surface of revolution with coordinates and velocity components. 
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Figurs 6. - Fluid particle with pressure and shear forces. 


COS 
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Figure 7. - Developed view of mean surface of 
revolution (fig. 5) in vicinity of R, 
showing relation between p , R d0 ^ and 
dR 


sin a 
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1.0 1.2 1.4 1.6 1.8 2.0 


R 

(a) Variation in with radius. 

Figure 8. - First group of numerical examples, showing effects 
of friction and heat transfer. 
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1.0 1.2 1.4 1.6 1.8 2.0 

R 


(b) Variation in pressure ratio with radius. 

Figure 0. - Continued. First group of numerical examples, 
showing effects of friction and heat transfer. 
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(c) Variation in flow direction with radius. 


Figure 8. - Continued, First group of numerical examples, 
showing effects of friction and heat transfer. 
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(d) Flow path for sin a = 1.0. 

Figure 8. - Continued. First group of numerical examples, showing effects of 
friction and heat transfer. 
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(e) Variation in polytropic, or small-stage, efficiency 
with radius. 


Figure 8. - Concluded. First group of numerical examples, 
showing effects of friction and heat transfer. 
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(a) Variation in with radius. 

Figure 9. - Second group of numerical examples, showing 
effect of diffuser wall spacing as affected by changes 
in compressor flow coefficient cp. 
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(b) Variation in static pressure ratio with radius. 

Figure 9. - Continued. Second group of numerical examples, 
showing effect of diffuser wall spacing as affected by 
changes in compressor flow coefficient q). 
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1,0 1.2 1.4 1.6 1.8 2.0 

B 

(c) Variation in flow direction with radius. 


Figure 9. - Continued. Second group of numerical examples, 
showing effect of diffuser wall spacing as affected by- 
changes in compressor flow coefficient cp. 


M 
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(d) Flow path for sin a = 1.0. 

Figure 9. - Continued. Second group of numerical examples, showing effect of diffuser wall spacing 
as affected by changes in coinpressor flow coefficient cp . 
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(e) Variation in polytropic efficiency with radius. 

Figure 9. - Continued. Second group of numerical examples, 
showing effect of diffuser wall spacing as affected by 
changes in compressor flow coefficient cp. 
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(f) Variation in vaneless diffuser efficiency (T]p )2 with compressor 
flow coefficient. 

Figure 9. - Concluded. Second group of numerical examples, showing 
effect of diffuser wall spacing as affected by changes in compressor 
flow coefficient cp . 
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(b) Flow path for sin a = 1.0. 
Figure 10. - Concluded. Design problem. 
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